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We present new results of our ongoing project toward a precision determination of the kaon B parameter with 
the Kogut-Susskind quark action in quenched QCD. New results taken at /3=6.4 and = 5.7 suggest that an 
apparently linear a dependence of Bk previously observed for /3 = 5.85 — 6.2 arises from a change of curvature 
from convex to concave as the lattice spacing is reduced. Fitting data for f3 > 5.93 with an 0(a 2 ) form yields 
Bk (NDR,2 GeV)=0.587(7)(17) in the continuum limit. We also describe a finite-size study of Bk at /3 = 6.0 and 
6.4, and a reanalysis of the theoretical argument for 0(a 2 ) behavior. 



1. Introduction 

Lattice QCD determination of the kaon B pa- 
rameter with the quenched Kogut-Susskind quark 
action has been pursued over a number of years. 
Last year we initiated a renewed effort toward this 
goal using VPP500/80 at KEK, improving upon 
a number of points of calculation of the previ- 
ous pioneering studies In particular, with 
a systematic scaling study over a wide range of 
lattice spacing carried out with improved statis- 
tics, we attempted to achieve a verification of the 
0(a 2 ) scaling violation theoretically suggested^ 
and then a precision determination of the contin- 
uum value of Bk- 



Our initial results taken at < 



= 1.3-2.6 GeV 



(/3 = 5.85 - 6.2), reported at Lattice 95 (J, how- 
ever, did not exhibit the expected 0(a 2 ) behavior 
but were apparently more consistent with an O(a) 
behavior. In order to understand this result, we 
have since extended our study by a new calcula- 
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tion both at a larger and a smaller value of a -1 . 
We have also carefully examined the question of 
finite-size effects in Bk, which could jeopardize 
the correct scaling behavior. In addition a re- 
examination of the theoretical argument for the 
0(a 2 ) behavior has been made. Here we present 
results of these analyses. 

In this article we mostly concentrate on data 
taken with an equal mass for d and s quarks. To- 
ward the end we briefly touch upon the question 
of quenched chiral logarithms that could affect 
the case of unequal quark masses. 

2. New runs and analysis procedures 

In Table [j] we list the parameters of all our 
runs carried out so far. Lattice sizes marked with 
* signify new runs since Lattice 95. Our main 
sequence of runs taken at a physical spatial lattice 
size of La r* 2.5fm now includes those at [3 = 
6.4(40 3 x 96) and at = 5.7(12 3 x 24). Finite- 
size studies are carried out at (3 — 6.4, and also 
at (3 = 6.0 to examine possible lattice spacing 
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Table 1 

Run parameters of our simulation. New runs since Lattice 95 are marked with *. 



a 


5.7 


5.85 


5.93 


6.0 




6.2 


6.4 




L 3 T 


12 3 24* 


16 3 32 


20 3 40 


24 3 64 18 3 64* 


32 3 48* 


32 3 64 


40 3 96* 32 3 96* 


48 3 96* 


#conf. 


150 


60 


50 


50 50 


40 


40 


40 40 


20 


skip 


1000 


2000 


2000 


2000 




5000 


5000 




er^GeV) 0.806(14) 


1.36(3) 


1.59(3) 


1.88(4) 




2.65(9) 


3.47(7) 




La (fm) 


2.9 


2.3 


2.5 


2.5 1.9 


3.4 


2.4 


2.3 1.8 


2.7 


m q a 


0.02-0.08 


0.01-0.04 


0.01-0.04 


0.02-0.04 




0.005-0.02 


0.005-0.02 




m s a/2 


0.0563 


0.0202 


0.0160 


0.0126 




0.0089 


0.0070 




fit 


6-16 


10-20 


12-26 


20-42 20-42 


14-32 


20-42 


25-69 25-69 


25-69 



dependence of finite-size effects. 

Gauge configurations are generated with the 
5-hit pseudoheatbath algorithm. Four values 
of quark masses, equally spaced in the interval 
given, are employed at each [3. 

The method of calculation of Bk is the same as 
reported at Lattice 950. We only summarize the 
main points here, referring to ref. [Q for details: 
(i) The physical scale of lattice spacing is set by 
the p meson mass in the VT channel calculated 
on the same set of configurations, (ii) For the case 
of an equal quark mass taken for d and s quark, 
we estimate half the strange quark mass m s a/2 
from mK/m p — 0.648. Interpolation of data for 
Bk to m s a/2 is made by fitting results for Bk 
to the form predicted by chiral perturbation the- 
ory, (iii) We employ tadpole-improved operators 
using the quartic root of plaquette as the mean 
link uq = P 1 / 4 . The one-loop renormalization 
factors, which include finite parts for the contin- 
uum MS scheme with the naive dimensional reg- 
ularization (NDR), are evaluated with g?L—(l/a) 
evolved from g^—^/a) = P/g% + 0.02461. (This 
point differs from the choice 9jfg( 7T / a ) adopted at 
the time of Lattice 950.) Setting \i = 2GeV in 
the logarithm of the renormalization factors, we 
extract Bk in the continuum MS scheme with 
NDR, which we denote as Bk (NDR, 2GeV). (iv) 
We estimate errors by a single elimination jack- 
knife procedure. 

3. Lattice spacing dependence of Bk 

In Fig. we present (NDR, 2 GeV) from our 
main sequence of runs as a function of m p a. The 
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Figure 1. Gauge non-invariant (circles) and in- 
variant (squares) Bjs:(NDR, 2GeV) as a function 
of m p a, together with a quadratic fit to four data 
points on the left (solid lines) and a linear fit to 
five points (dashed lines). 

left and rightmost points are the new data taken 
at P = 6.4 on a 40 3 x 96 lattice and at @ = 5.7 
on a 12 3 x 24 lattice. 

As we observed at the time of Lattice 95, the 
middle four points are consistent with a linear 
a dependence. Viewed more closely, the fourth 
point on the right lies below, while the other three 
are on a straight line. This suggests a convex 
shape of the curve of Bk toward larger lattice 
spacings, which is confirmed by the new data at 
(i = 5.7 (the rightmost point in Fig. |l|). 

On the other hand, the new data at j3 = 6.4 
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Figure 2. B K (NDR, 2GeV) as a function of phys- 
ical spatial size La (fm) at (3 ~ 6.0 and 6.4. 

at the leftmost lies higher than the straight line, 
which is suggestive of an onset of an 0(a 2 ) be- 
havior for smaller lattice spacings. The solid line 
shows a fit of form Bk = cq + ci{m p a) 2 to the 
four points, excluding the two points on the right 
at (3 = 5.85 and 5.7. 

We note, however, that the data by themselves 
still do not allow a distinction between an 0(a) 
and 0(a 2 ) dependence; five points for m p a < 0.6 
can also be fitted with the form Bk = cq + cim p a 
(dashed lines). 

Let us also note that the discrepancy of values 
between gauge invariant and non-invariant oper- 
ators, even in the continuum limit, represents a 
systematic uncertainty of 0(g i ) arising from the 
renormalization factors which are correctly incor- 
porated only to the order of one loop. 

4. Lattice size dependence of Bk 

We should emphasize that a distinction be- 
tween an 0(a) or 0(a 2 ) dependence rests on a 
very tiny effect in the range of lattice spacing 
being examined; a 2% decrease of the point at 
[3 = 6.4 would render the 0(a) behavior strongly 
favored. Clearly a reexamination of finite-size ef- 
fects is needed, especially since our physical spa- 
tial size decreases by about 10% from La = 2.5fm 



to 2.3fm between f3 = 6.0 and 6.4. 

In Fig. ^ we summarize results of our finite- 
size study at f3 = 6.0 and 6.4, plotting Bk (NDR, 
2GeV) as a function of physical spatial size La. 
For La > 2.3fm, we do not observe a size- 
dependent shift of Bk beyond the statistical error 
of at most 1%. Therefore, we conclude that the 
data in Fig. [|, taken with La=2.3-2.9 fm, are not 
affected by finite size effects. In particular, the 
rise of Bk at (3 = 6.4 is not a finite size effect. 

An unexpected feature in Fig. || is that finite- 
size effects below La ss 2fm is markedly enhanced 
at a weaker coupling of (3 = 6.4 than at 6.0. Ori- 
gin of this behavior is not clear to us. 

5. Reexamination of the theoretical argu- 
ment for 0(a 2 ) behavior 

A subtle point in the theoretical analysis of a 
dependence of Bk is that numerical simulations 
employ 4-component hypercubic fields q(y),q(y) 
defined on a coarse lattice of lattice spacing 2a, in 
terms of which the Kogut-Susskind quark action 
has an apparent 0(a) term: 



S 



KS 



V 



q(y)[(l P ®I)D p + jn(I®I) 



+a(75 ® ^5)D; + aigF^T^/A 
+0(a 2 )]q(y). 

However, a redefinition of fields given by 

Q = (i - p M 2 )i> Q = «(! - E p » / 2 ) 

with = (1 - 7^75 (g> £m£5)/2 removes the 0(a) 
terms at tree-level so that the action becomes 

V 



+0(a 2 ). 

For bilinear and four-quark operators it is also 
possible to show that 

(P\0(QQ)\P') 



[1 



i(P' p -P p )/A](P\0(qq)\P') 



+0(a 2 ). 

Therefore Bk evaluated with q(y)'s with P = P' 
has no O(a) corrections at tree level. 
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The rest of the argument essentially follow 
that of ref. j|. In the presence of interac- 
tions the action expressed in terms of Q fields 
may contain 0(a) terms. However, using sym- 
metry of the original action under translation, 
Q(V) {I®^)(l+dn/2)Q(y'), rotation, Q(y) 
((1 - 7m7«0 ® (1 - M)/2Q(y'), and U(1) A 
symmetry, Q(y) — » e la ^ 5 ^^Q(y) performed 

simultaneously with a mass rotation, m —> 
e -m(75®?5) TOe -M*(75®? 5 ) ; one can show that there 

are no dimension five operators giving rise to an 
0(a) term. Therefore, the original action is itself 
improved up to 0(a 2 ). One can also check that 
the Bk operator is similarly improved to 0(a 2 ) 
since there are no dimension seven operators with 
the same symmetry properties as the original one. 

6. Extrapolation to the continuum limit 

Our theoretical reexamination supports the 
original assertion|| that there should be no 0(a) 
corrections in Bk including gauge interactions. 
This leads us to attempt a continuum extrapola- 
tion of our data assuming the 0(a 2 ) dependence. 

Fitting data in Fig. |l| to the purely quadratic 
form Bk = cq + ci(m p o) 2 , excluding the two 
points on the right at (5 — 5.85 and 5.7, we find 

Bff (NDR, 2GeV) = 0.587 ± 0.007 ± 0.017 

where the first error is statistical and the second is 
the systematic error of 0(g i ) estimated from the 
difference of extrapolated values for gauge non- 
invariant and invariant operators. Changing the 
number of points and/or including higher order 
terms in m p a, we find that the extrapolated val- 
ues stay consistent with the above value. 

7. Bk with unequal quark masses 

We measure Bk also for unequal quark masses 
for d and s quark in our simulation. In this case, 
a potential problem of quenched chiral logarithm 
may arise in the extrapolation to the chiral limit 
in m d . 

We examine this problem through fits of our 
data for rrid ^ m s to the form predicted by 
quenched chiral perturbation theory Gener- 
ally we find very good fits. However, the param- 
eter S = TOo/487r 2 / 2 multiplying the quenched 



chiral logarithms turns out negative albeit with a 
large error, 5 « —0.3 ± 0.3 (we note that the sign 
of the 5 term in ref. || is incorrect || ) . Also fits 
excluding the logarithm are quite acceptable, only 
slightly worse than those with it. Thus we cannot 
conclusively conclude the presence of quenched 
chiral logarithms. 

Let us add the remark that the mass inside the 
quenched chiral logarithm should theoretically be 
that of the non-Nambu-Goldstone pion||. Fits 
in this case differ from those with the Nambu- 
Goldstone pion mass only very close to m d = 0, 
however, so that the conclusion above does not 
change. 

8. Conclusions 

Our new data for Bk extending the results of 
last year to smaller and larger lattice spacings 
suggest that the curve of Bk changes curvature 
from convex to concave as the lattice spacing is 
reduced from m p a w 0.9 to 0.2. This behavior 
explains an apparently linear dependence we ob- 
served in our previous data, and is also consistent 
with an 0(a 2 ) dependence theoretically expected 
toward the continuum limit. To confirm if Bk 
stays on an 0(a 2 ) curve, we plan to carry out one 
more run on a 56 3 x 96 lattice at (3=6.65 with yet 
smaller lattice spacing of m p a w 0.16. 
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